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BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 5, 1943 


OUTLINE OF A NEW MATHEMATICAL APPROACH 
TO GENERAL BIOLOGY: II 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


In this continuation of a previous paper another formal principle 
is added, similar in its form to Hamilton’s principle in mechanics, and 
based on a mathematical formulation of a postulate, stated by A. J. 
Lotka, namely that the events in the living world are determined by the 
requirement of the maximum of the total flow of energy. 

The principle is applied to a wide number of biological phenomena, 
including paleontologic developments, relations between death and birth 
rates of organisms, ontology, cell division, ete. 


II. The Organic World as a Whole 

The environment of each organism being principally composed of 
other organisms, the organic world as a whole is characterized by an 
incessant interaction between different organisms. Some of them 
provide food for the others. Some are detrimental to others either 
by feeding on them or inhibiting the activities of these others by 
metabolic by-products (bacterial toxins). Thus the whole organic 
world is a system in which continuous exchanges of matter and energy 
takes place, a system of complex flows of matter and energy. 

A. J. Lotka (1922a, b) suggested as a possible general biological 
principle, that the processes of life run in such a way as to maximize 
the total flow of energy in the organic world. We shall attempt to de- 
velop this suggestion further, and formulate it in mathematical form. 
We shall first proceed in a purely abstract fashion. 

Let any organism be completely determined by n determining 
parameters 

U1» He, H3°°* Un. (1) 


Consider first the case that those parameters have only a finite set of 


discrete values 
te , vo , g(t «++ ay 


ads ies hae Vy 
(2) 
em .. (m) : 2m) : 16) 
where in general m >, and ™ is in no way related to n .. Subsequent- 
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ly we shall pass to the continuous case, letting all x,’s vary continu- 


ously. 
Consider an organism, that is characterized by the following set 


of determining parameters. 
a é A ; x8 +5 Uhm (3) 


and let the total number of such organisms in the organic world be 


Keiyieaeactoe es (4) 
Let another organism be characterized by a set of values 
x) ’ Lo"? , x3) po ,t* ’ (5) 
and the total number of organisms of that type be 
ee a. (6) 


We may denote the energy exchange of interaction of the first 
kind of organisms with the second kind by 


DLR tonks iN pers Ne (7) 


The coefficient a, is to be considered as positive if the interac- 


12K} U4...Un 


tion is such, that it increases the energy content of the first type of 
organisms, negative—otherwise. 
The total energy of interaction or total exchange in the whole or- 

ganic world is then given by the sum 
Loe Se Cee ey es N ° (3) 


ki; I 14...0n 


We may also consider the interaction of each type of organism 
with the inorganic world; this will-be of the form 


ee > sare Ne c (9) 


It may be remarked, that we do not need necessarily to have 
Qin = — Mi, aS One might be inclined to assume at first. The anti- 
symmetry would be required only if the organisms would interact in 
isolated pairs. In general, when there is simultaneous exchange of 
energy between many organisms, a3, = — dy; . 

We may now consider different types of postulates that lead us 
to interesting equations, governing the development and behavior of 
the organic world. 

First we may postulate that the configuration of the whole or- 
ganic world tends to such a configuration which maximizes the value 
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I, + I,. When such a configuration is reached, no further changes 


occur, and if F UN a) ig any function of Ne ; 


dl (N,. 3.) UNG es 
Saar ie ana (Ve) eis cae 0! (10) 


where F” denotes the derivative. 
At the same time, when J, + J; has a maximum 


@ (Tents) 


CINE x Cae 
Hence, it is natural to put 
dF 
(Ni sen) ao sali) (12) 
dt a CNiere ‘ 


where F is to be determined later by biological considerations. Intro- 
ducing equations (8) and (9) into (12) we find 


dE (N,. in ) 
ter dipuna Oe Nes ze obs 2 (13) 


If we now let the values of x; vary continuously, the coefficients 


a eee become functions of the values 21; 2 --- % of the first or- 

ganism and of the values x,', x.’ --- x,’ of the second. The coefficient 

Rs , become functions, of 41, X +--+ Xn, N.,, ,, becomes a function 
a 


N (41+*+ &n) Of 41; Lo+++ Xn, Nos becomes a function N (a1 +++ Xn’) of 
the x;’s and the system of linear differential equations (13) becomes 


an integrodifferential equation of the form 


ON (1, +++ Gn, t) 
ot 


= S fives SK (Gq, +++ Ln 5 Hr! +++ En’) Ney Ey OR a, (14) 
+ f (%1, %2°-+ In), 


the integration being extended over the whole range of variation of 
the #;”s. 

If the kernel K (2, , 2 +++ Unj 1’ «++ %n') and the function f(x, , x2 
-++ %n) are known, the integrodifferential equation (14) determines 
for given initial conditions, N(a1; «++ tn, €) = No( a1, ++ &n) for 
t = 0, the function N(x,,--- x, ¢) for any time t > 0. Since the 
values of 2, , 22 -:: %, determine the shapes and sizes, and in general 


F'{N (415 +++ &n, t) } 
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the anatomic structure of the organism, the knowledge of N(%1, %2 
.»+ 2, t) gives us thus the distribution of different types of organisms 
in the organic world at any time t > 0. For certain values of the x4’S 
the function N(x,,--- x,, t) may be zero. Organisms with anatomic 
structure corresponding to such values of the determining parameters 
will not exist. As we have remarked elsewhere (Rashevsky, 1941), 
the knowledge of the distribution function N(2#,,--- %,, t) gives us 
also the correlation coefficient between the different determining para- 
meters, if nm = 2. For practical purposes, whatever we may wish to 
obtain from equation (14), it will be rather difficult to do for such a 
large number of determining parameters as is likely to actually occur. 
We may however fix our attention on any pair x,, x, of the n deter- 
mining parameters, and wish to find the correlation between them, 
when the remaining » — 2 determining parameters vary at random. 
To this end we may integrate the function N (#1, %2 ++: %,) with re- 
spect to all the variables except x, and x;, the integration to be ex- 
tended over the whole range of variation of those n—2 variables. Thus 
we obtain 


N (Gia, %,),) = f---fN (Hr +++, Un, t) dx dte 


oo dda Db par ie d2l,4 Ofee rae! ce ° eA. 


Since, denoting by N(t) the total number of organisms in the world 
at the time t , we have 
Sos SN (a1 +++ Xn) Ax, +++ dt, = N(t) , (16) 
therefore 
SIN (a, , &,, t) dx, dx,=N(E). (17) 
Hence N(x,, «,;, t) may be considered as a distribution function for 
+, and x, for random values of the other determining parameters. 
By introducing similarly a function 
Fi (Ooo diy wee) 
= So fK (ay se My Uy +88 Cy OU -* OF see 
ADs Aer +++ Ty +++ dity' +++ Me yg OM pay ++ OH py AL guy ++ AL'y 
and by applying the same reasoning to N(x,, x,) in regard to mutual 
interaction, we find 
ON (a,, He,t 
F'{N (a,, %s, t)} Aaa th 


= [SK (trom.) Nae, &)) dada Flap, a, as 
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where 
Preps ie hf (tyee8 ae) da, + 


Ady Ars Bis AXe4 1 teak Oa;, . 

Equation (18) may be used for the determination of the coeffici- 
ent of correlation between x, and x, and to the extent that those may 
determine certain features of the anatomic structure, it may be used 
for a theoretical prediction of the correlation between some anatomic 
structures. : 

Equation (14) describes the gradual development of the whole 
organic world. The change of N(a,,--- x,, t) in time gives us the 
change of the general structure of the organic world with changes 
from one type of organism to another, as is found in paleontology. 
One of the problems of the theory is to determine K (2, +++ &nj.X1'+++ Hn) 
and f(x, --- x,) from theoretical considerations, and to see whether 
this will give us such an N (a, --- 2, , t) which would roughly describe 
paleontological facts. 

If f(%, +++ &,) > 0 for any values of x, -:: %,, equation (14) im- 
plies that even when N(x, --- %,) = 0, still for every set of values of 
U1 +++ in, ON (Xi +--+ Xn) /ot > O. This would mean, that organisms of 
any type could be spontaneously generated from nonliving material. 
Since this is not the case f(%1, --: %,) cannot be positive everywhere. 
Neither can it be negative or zero everywhere, since then the organic 
world could never have started, because if f(%, --- %,) = 0, then for 
N (#1, ++: tn) = 0, ON (a, --: &,) /ot = 0. At some geological epoch 
life must have started out of the nonliving. And it did start with the 
simplest organisms for which most of the x;’s are zero, and the others, 
perhaps only M and q, rather small. To obtain this result we may 
assume f (21 --: %,) to be of the followirfg form. 

Let 

%,—=M,2.=—¢q, and let A, e, and « be small positive quantities. 
Then let 

f(y sa <> fp) — A for 0< 2, < &,,0 <a. < &, te = Ue 
— 0 

Ge tend, J — 0 fOr any; Ot in (19) 

With such a form of f , we shall find that even when N (2, --- %n) 
= 0, dN (a, --- %,) /at > 0 for such organisms, which correspond to 
0<4,<a4;0< 4, <a. But dN (2, --- x,) /dt = 0 for all higher or- 
ganisms, which have any of the x; (¢ > 2) greater than zero. 

Such a situation would still imply, that even at the present time 
the spontaneous generation of some lowest forms of organisms takes 
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place. Though there is no direct evidence against that, if such an evi- 
dence were forthcoming, it could be met in two different ways. 


A. The first and simplest way is to make f depend on one or more 
parameters 7(t) which vary with time. Those parameters may char- 
acterize the physico-chemical properties of the environment, such as 
temperature, pressure, humidity, etc. at different times. We may now 
put instead of (19), a different requirement, namely: 

Requirement (19) holds only for 4 < mo; for > mo f(4%, %2 
-++ %,) = 0 for any values of «;. If now 7 increases with time, being 
originally less than 7, , then after a while spontaneous generation will 
become impossible, for any kind of organisms, no matter how simple. 

B. The second way is to make f not-a function of the x;’s , but 
a functional of N(x, --- %,) ; in Volterra’s notations f|[N (a1 --- %n)]|, 
such that 


F\LN (a1 +++ %n)]| > 0 


when N(x, -:- %,) is different from zero only for 0 < 4% < 4&3 
0< %. < fos le — ee 


flEN (5+. 25) ])/=0 


for any other shapes of N(2,, --- %n). In that case, f will become 
zero aS soon as the organic world develops sufficiently to have even a 
few organisms with x; > 0 (i > 2). 

There is, however, a good reason why the proposed equation 
should in principle not be adequate to take care of the problem of the 
origin of life. This reason we shall discuss in a subsequent paper. For 
the time being we shall keep the f term, since an exchange of energy 
with environment always takes place, and take care of the general 
non-occurrence of spontaneous generation at the present time by a 
proper choice of the function F’. 


The equation 
K (1 5 +++ hn; Hy’ +++ Ln’) =0 (20) 


defines for a given set of values x, --- x, , considering x;’s as variables, 
an  — dimensional hypersurface, which divides the whole space of 
the variables x;’ in two regions, one for which K > 0, and another 
for which K < 0. The hypersurface (20) may be multiply connected. 
Let us denote by V. the region in the hyperspace for which K > 0 4 
and by V_ the region for which K < 0. Since N(a, --- 2, t)= 0, 
therefore we may write equation (14), thus, denoting by vertical 
lines absolute values. 
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ON (#1 +++ Un, b) 
ot 


Seer S|K (a. ++: Ln> Ly +++ Xn) |N (a1 
"++ Vn) Eats’ +++ Ein! ase JESS iL KG (rye arn ar 
Nigie. a, dade. +f. 


both integrals being positive. 

Considering 2, , 2 --- %, as fixed, we may interpret the first inte- 
gral of (21) ag the increase of organisms of type (%1, 22 --: %n) by 
birth, while the second integral may be interpreted as the decrease of 
the number due to death. Both birth and death rates thus appear as 
depending on the presence of other organisms. For such values of the 
x;’s, for which f > 0, the function f adds to the birth rate, while for 
such values for which f < 0 it adds to the death rate. To bring equa- 
tion (21) into the usual form, when the birth and death rates are pro- 
portional to the population, we must put 


F'{N (21, +++ Ln) } 


(21) 


aes 2°) 


F{N (41 +++ %n) } =log N (41 -++ 4n). (22) 
Equation (21) then becomes 


ON (x1 ote Uns t) 
ot 


= N (2; *-- an, €) Syoet S| K (ats ++ ttn5 tr! + ata’) |N (ats! 
228, ) day =~ br," 

SN Gao ase fie (arn: 0 tn) |N (aa 
Bae Ly) OL Oy ON | 


(23) 


This may be considered as a more general form of the usual 
growth equation: 


dn ‘ 
ia an — bn 
which may be written in the form 
‘ oan — (b' + c'n)n 


with 
a=—a—b'; b=c’. 


Thus if the kernel K is known, by means of equation (20) we may 
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determine the birth and death rates for any species characterized by 
a given set of values x, --- 2, of the determining parameters. The 
knowledge of those quantities gives us also the average lifespan of 
the species. If the death rate is equal to g(t) and m denotes the total 
number of individuals at a given time, then, in the absence of birth 


dn 

hath 2s 24 

dt q(t)n (24) 
or 

N= Noefut ai (25) 


which shows that ” tends to zero exponentially. 
Let the distribution of life-spans of a group or organisms be such, 
that the number of individuals that die at the age t, is given by 


MoN (t) at 
with 
[ X@)at= iv (26) 
We have for the average f 
z= | “tnctydt. (27) 
0 


If we consider that originally all individuals have the same average 
age, then in absence of births, the number of individuals, originally 
N> , becomes at the time f: 


ey [ x@yat=n [seat (28) 
t 0 


Let at the time ¢*, this number drop to 1/e of its original value. Then 
we have 


if "Neat =} (29) 


e 
which solved for ¢* gives, because of (27) 
*=u(t), or £=“a(t*). (30) 


From (25) we have for the time when n, drops te 1/e of its origi- 
nal value ; 


{vaat= 1 (31) 


which determines ¢* for a known q(t). Substituting that value of t* 
into the second equation (30), we find é. 
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If the distribution of the life spans is due to chance, N(t) will 
be the normal distribution and ~ may be found explicitly in terms of 
the parameters of that normal distribution. 

Thus, knowledge of K in equation (23) gives us the birth and 
death rates and the average life span for any organism. 

It may appear strange that both birth and death rates are depen- 
dent on the presence of other organisms. Concerning the birth rate 
this is however quite natural. For the multiplication of organisms 
goes on only under proper conditions of nutrition, which involve the 
destruction of other organisms. Even as regards some simplest micro- 
organisms, which feed on non-living media, they do not multiply in a 
culture where only a few organisms are present. 

In regard to death rate the conclusion is no less natural. Equation 
(23) gives the actual death rate, which is due to all possible causes. 
And with rare exceptions, death occurs due to interaction with or- 
ganisms. Death of pure old age is exceedingly rare, even amongst 
humans. If f > 0, then in terms of our equations (23) we may ex- 
press that death rate by taking the integral 


[tees fe (igen ne Pies 0, ') da,’ ++> dx, (32) 


over a very narrow strip of values in the immediate neighborhood of 
the hypersurface (20). If 7, is a small quantity, we may integrate 
over the range 


Gee Fee U/l ’ ve sf 1 ’ (33) 


x;° being a coordinate of a point of the hypersurface (20), and 7, be- 
ing considered as the range of biological variations of the character- 
istics of the species. Expression (32) integrated over the range (33) 
may be divided again into a positive and negative part, the latter giv- 
ing us the death rate from purely old age. Due to the smallness of 77; , 
it will be very small. 

The problem of handling equation (23) resolves itself into two 
parts: first the problem of determining from physico-chemical and 
biological considerations the form of K (a1 , +++ %n; X1' +++ Xn’); and sec- 
ond—the solution of the nonlinear integrodifferential equation (23). 
The similarity of equation (23) to some nonlinear systems of differ- 
ential equations with many variables studied by V. Volterra (1931) 
suggests the approach to the second problem by a transition to the 
continuous case from Volterra’s equations. The first problem opens a 
very wide field for further theoretical studies. 


It is important to remark, that the limits of integration in equa- 
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tions (14) and (23) are not all constant. As we have seen in part I, 
if M and q are prescribed, the velocity v may have in general an up- 
per limit, which will be a function of M and q. Thus if we deal for 
instance with only three determining parameters 7%, = M, «2, = q and 
23 = Vv, we may integrate with respect to x, and x. from 0 to o, but 
with respect to x; only from 0 to a value v(M,q). In other words the 
total space is limited by the planes g= 0, M=0, v = 0 and the sur- 
face v =v (M, q). 

Limiting ourselves, for simplicity, to only two determining para- 
meters, M and n, we may for instance determine the form of K by 
such considerations. For any small value of 1 the interaction between 
the two organisms M and M’ will be such that the greater the ratio 
M/M'’, the greater the interaction. (Remembering that by definition the 
interaction of M and M’ is positive, when it favors M.) This will hold 
however only within a certain range of ratios M/M’. A too large ra- 
tio M/M' may be unfavorable to the first organism. 

Thus while an elephant may readily have the better part of a 
tiger, yet it may be quite helpless against some pathogenic bacteria 
with which its interaction will be negative. Hence K must for a fixed 
n vary in such a way as to have a positive maximum for a certain 
M/M’, then drop again and become negative. 

But as 7” increases, the advantage of a larger mass becomes less 
and less important. A smaller animal with greater nervous coordi- 
nation may have the better of a larger one with a less developed ner- 
vous system. Denoting by A , B, C and D four constants, we may try 
the following expression for K(M,n; M', 7’): 

(A +Bn-—) — (A + Bi 1) 

K(M,n; M',n’) Con 

Da 


D+ 


(34) 


Expression (34) has the property, that for fixed M and n and for 
very small values of n, the first term is the most important. With in- 
creasing n , not only does the first term decrease relatively to the sec- 
ond, but the optimal value of M/M’ shifts towards larger values of 


M/M’. For very large n , the second term becomes the only important 
one. 


Without solving equation (28), we may surmise what conse- 
quences of such a choice of K will have. Given an initial condition, 
when only small values of M and n are present, we shall find that the 
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larger the value of M, the larger the birth rate and the smaller the 
death rate. Hence dN (M , 1, t) /dt will be the largest for large values 
of M , and the value of N(M , n) will become large for large M’s. As 
however the value of N(M, n) increases for all values, and therefore 
the number N(M, 7) of organisms with large n also increases, the 
birth rate will become largest for organisms with a large». Such or- 
ganisms begin to prevail and those with large M , but small n actually 
die out. (C.f. preponderance of large reptiles in Jurassic era and their 
recession. ) 

As a final generalization of our approach we may suggest, that 
instead of postulating the maximizing of the expression J, + I,, as 
given by expression (8) and (9), we may introduce a Hamiltonian 
principle and postulate that (Whittaker, 1927, p. 249) 


t 
J ig - $f D> Dem é Nt) =: (35) 
ta 
With 
1 AN. on 
Beate ls Biota Nos dine (36) 


requirement (35) leads us to a system of non linear equations, which, 
in passing to the continuous case lead again to equation (23). 


TI. Further Applications; General Discussion 

Equation (23) or the more general equation (14), represents a 
theory of the development of the whole organic world. Known facts 
of paleontology should be deducible from those equations in principle, 
provided we can overcome the mathematical difficulties of solving 
them. ; 
We may apply the same principle, which leads to equation (14) 
or (23) to a number of more restricted, though no less important 
problems, by generalizing it somewhat. We may assume that the 
principle (35), or the requirement of the maximum of I + I,, holds 
not only for the organic world as a whole, but for any part of it. We 
may apply it then either to limited organic groups, or even to indi- 
vidual organisms or cells. 

An example of such a limited group is given by human society. 
While the gross physiological properties of man may be considered 
as rather uniform constants, different psychological traits vary con- 
siderably from person to person. The social interactions of two indi- 
viduals are determined by the values of the parameters which meas- 
ure those psychological traits. By exactly the same argument as in 
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section IJ we shall arrive at an integrodifferential equation of the 
form (23) for the determination of the distribution function of dif- 
ferent traits. Not only does this equation give the distribution func- 
tion, but it also gives its variation with time. As has been shown 
elsewhere (Rashevsky, 1942) such a distribution function determines 
the social structure of society, and its variation with time. Thus 
mathematical sociology actually becomes a branch of mathematical 
biology, and history may in a sense be considered as a small branch 
of paleontology. 

The distribution of different psychological types may also be 
studied in this way. The methods of theoretically predicting the cor- 
relation between any pair of determining parameters, mentioned in 
section IJ, suggests a possible approach to a theory of the experi- 
mentally found correlations between Thurstone’s primary factors. 

The same procedure should lead to a theory of distribution of 
different pathological variations, giving us the incidence of certain 
organic diseases. 

‘All human activities may be considered as biological manifesta- 
tions and the principle of maximum energy flow may be applied to 
all of them. If so, then even such problems as the distribution of 
wealth may be treated in a similar way. Wealth is an important fac- 
tor producing energy flow, and the Pareto distribution may perhaps 
be derived from considerations of maximizing that flow. 

A characteristic feature of equation (14) or (23) is that they 
imply the necessity of death in general for all organisms, if the or- 
ganic world is to function as a whole with a maximum flow of energy. 
This “necessity of death” may be considered also from another angle, 
in reference to humans alone. Destruction of one’s own fellow men 
is a process that involves itself a large expenditure of energy and 
hence produces an energy flow. It may be, that actually the energy 
flows involved in such a destruction or war, are greater per unit time 
than those involved in peace time activities. A war however cannot 
last indefinitely, because it would lead to the extinction of the race. 
But a maximum energy exchange per unit time may be obtained by 
intermittent alterations of war and peace. 

For simplicity let those alterations be assumed quite regular. 
Let the duration of peace be t,, that of war t,. The rate of energy 
exchange during peace is an increasing function of time since the 
population increases. The rate of energy exchange during war is a de- 
creasing function of time, since gradually the population decreases. 
The total change of energy during t, is f,(t,), the change during ft, is 
fw(t»). The average exchange per unit time is 
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—_ Folte) + fro (tr) 


Oty to + ty (37) 


Depending on the form of f, and f», fa may have a maximum for 
a set of values t, , t,.. Those values determine the corresponding aver- 
age duration of peace and of war, while t,,/(t) + t,) determines the 
frequency of incidence of wars. 

In passing to the considerations of individual organisms, an in- 
teresting possibility at once suggests itself in regard to general inter- 
nal structure and ontogenetic development of the latter. 

The different organs of an organism interact with each other, the 
interaction resulting in energy flow towards or away from the organ. 
Let us denote by M, the total mass of the organism, by M; the masses 
of different organs. We may consider as a first approximation the 
energy flow from organ & to organ 7 as given by 


dir M; M;, . (38) 
Moreover the interaction of the organ i with the whole organism may 
be put in the form 
Qin M;, My. (39) 
The interaction of an organ with the external environment may be 
taken as 
b, M;. (40) 


The growth of an organ results in an exchange of energy which in 
general will be of the form. 


c(M;) ae (41) 


We may now introduce a Hamiltonian principle for the organism, 
of the form 


[OF +E Q0Mo =) 


with 
F=>>4%MiM,+>60:M;; 
4k i 


1 dM; 
a ee iM, 


This will give us a system of differential equations, determining 
the variation of different M;’s. As initial condition we may put all 
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M; zero except M,. This is the condition of an undifferentiated ovum. 
The gradual development of the different organs should then be de- 
scribed by the system of differential equations. 

By a similar argument as in section IJ, we may limit our consid- 
erations to only a few more important organs, such as alimentary 
tract, circulating system, respiratory system, muscles, and nerves. 
The values of a’s may then either be calculated from theoretical con- 
siderations, or determined by comparison of the integral curves of the 
system of differential equations with actual embryological observa- 
tions. 

The variational principle may be applied also to problems which 
have already been dealt with in mathematical biophysics. In the the- 
ory of the central nervous system we usually introduce hypothetical 
structures, consisting of excitatory and inhibitory fibers, and invent 
more or less those structures ad hoc to explain certain phenomena. 
We may however proceed now in a more satisfactory way. Develop- 
ing first an abstract theory of neuronic network, along the lines initi- 
ated by A. S. Householder and Walter Pitts, we may then put the 
problem in this way: For a given sequence and type of stimuli, which 
of the theoretically possible structures will give rise to the greatest 
energy interchange? A structure adopted on this principle will be 
free from the ad hoc nature, inherent in structures studied hitherto. 
It may perhaps be generally proven, that any such structure will have 
to contain both excitatory and inhibitory fibers, and thus the intro- 
duction of the latter will cease to be an independent hypothesis. 


If excitation results in energy release, then our concept of pleas- 
ure as a general excitatory phenomenon and pain as a general inhibi- 
tory phenomenon may be reformulated thus: any situation which in- 
creases the total positive energy flow within the organism is pleasant, 
any situation that decreases it is unpleasant. 

Application of the extremum principle to individual cells may 
be illustrated on the problem of cell division. 

Using H. D. Landahl’s (1942) notations, we shall find in general 
that the energy flow, both inside the cell, and between the cell and tts 
surroundings, is a function F(e; 7, = 5 “) of e, r, de/dt and dr/dt. 
The quantity de/dt = e¢ enters into the function F in two different 
ways. First of all, during the process of elongation, as the surface of 
the cell increases, there is a transformation of some other form of 
energy into the additional surface energy. The rate of that transfor- 


mation depends on de/dt = «. Second, the dissipation of energy due 
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to viscous friction depends on <. In a similar way dr/dt = r is intro- 
duced. 


By requiring 
te 
6 Fdt==0 
ti 
we find the system of simultaneous differential equations 


a (OF _ oF 
dt de de 

42 
@ ak, ar Ses 
dt or or 


We thus describe the elongation and constriction by one set of 
equations, instead of treating them separately, as hitherto. This must 
be considered as a definite advantage. 

The energy flows include such things as the dissipation of energy 
in the diffusion fields, as well as the energy transformations in various 
metabolic reactions. The calculation of the former will involve the 
diffusion coefficients, permeabilities, and reaction rates. The calcu- 
lation of the latter will involve again the reaction rates and also the 
heats of reactions. Thus the cell division curves will depend also on 
the latter, as one naturally should expect. The missing of the heats 
of reactions, and of the general energetics of the cell, from the pres- 
ent theory of cell division has been one of its defects. The whole prob- 
lem is treated in a separate paper. 

As a last illustration we may apply our general principle to the 
problem of frequency distribution of words of different lengths. A 
spoken word results in certain activities, overt or covert, and thus 
speech produces energy transformation. Let the number of words of 
length 1 be N;.. The energy transformation produced by a word will be 
in general a function f(/) of its length, that for purely physiological 
reasons. But the energy flow created by a spoken word will also depend 
on the frequency of use of that word in the speech, and be a function 
U(N,). In the collection of words of different lengths, the total energy 
flow will be 


l=00 
> Nil[f(l) + U(N:i)]). (43) 
1=0 

The requirement of maximum of (43) gives 


dU 
F(t). -.U (Ni) EMG Ta 0. (44) 
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If we assume as an example only 
f(l) =alogl, 
U(N,) = blog N:, 


which may be made plausible by psychological considerations, equa- 
tion (44) gives 


(45) 


alog!+ blog Ni+b=0. (46) 


G. Zipff (1935, p. 23 ff.) has studied quantitatively the relations be- 
tween the frequency of occurrence of words and their lengths, in dif- 
ferent languages. Such data may be used to. verify such or similar 
relations as the above. 

The author is indebted to Dr. Alston S. Householder and to Dr. 
Herbert D. Landahl for reading the manuscript and valuable com- 
ments. 
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A suggestion is made to establish a general variational principle 
of which the recently discussed principle of maximum energy exchange 
in biology and the usual Hamiltonian principle in physics would be par- 
ticular cases. 


In a previous publication (Rashevsky, 1943), we have formulated 
mathematically a hypothesis stated some time ago by Alfred J. Lotka 
(1922 a, b), namely, that the course of events in the biological world 
is determined by maximizing the total energy exchange between differ- 
ent biological units. We have suggested that this principle may be 
stated in form of a variational principle. If x;’s are any kind of bio- 
logical quantities that determine the energy exchanges in a, biological 
system, and w,’s are the corresponding derivatives with respect to 
time, we construct a “Hamiltonian” F'(a,,--- %,,.@1 +++ 2), and re- 
quire S 


i ae 
6 | Fat=o, (1) 


so as to maximize the integral. Numerous possible applications of 
that principle to a great variety of biological problems have been 
discussed in loc. cit. 

In physics, the “Hamiltonian” H is a function of the generalized 
coordinates gq and momenta p , and we have 


s | Hat=0, : - (2) 


minimizing the integral. 

The question naturally arises as to whether principle (1) may 
not be reduced to (2), so as to avoid a duality of living and non-living. 
This does not seem likely, and A. J. Lotka (1922 a, b) has proposed 
his hypothesis as an independent principle. It must be remarked that 
F has the dimension of energy flow, hence [F] = mi? t*,, while H has 
the dimension of energy, and therefore [H] = ml? t?. One might per- 
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haps think that in some way both principles operate simultaneously, 
both in the living and non-living, but that the role of the first prin- 
ciple is more pronounced in the living world, that of the second in the 
non-living. The most natural suggestion would be to put 


be 
of (a, F — a, H)dt=0, (3) 
ty 


and assume that the first term prevails in the living, the second in the 
non-living. If, however, we make such an assumption, we must as- 
sume a, and o to be of different physical dimensionality, so as to make 
a,F' and a.H of the same dimensionality. This introduces a rather un- 
pleasant indefiniteness. 

It may be suggested that a, and a, be chosen so as to make the 
whole integral in (3) a pure number. This would give 


Roa aay ee ih AO rel PS ee (4) 


We may make the hypothesis that a, and a are two universal 
constants. All systems for which o,f >> a.H we may define as bio- 
logical; all systems for which o,f << a,H we may define as physical. 
For each of those classes both a, and a, drop out of the picture. But 
when oF is comparable to a.H , we have intermediate systems in 
which a, and a, play a role. 

However, inasmuch as mathematically only the ratio 7 = a,/o2 is 
important in the variational problem given by equation (3), therefore 


it may be better to leave the integral of the dimension energy X time 
by putting 


6 { G@r-mat=o. (5) 


We may assume 7, which has the dimension of time, to be a uni- 
versal constant characteristic of biological phenomena. If 7 is very 
small, then systems in which energy exchanges are not too large fall 
into the domain of physics. Systems in which the energy exchanges 
play a preponderant role fall in the domain of biology. When 7 F is 
comparable to H_ we have borderline systems. Such may be perhaps 
the case of single cells. By developing a theory of such systems, we 
may determine the constant 7 from any comparison of quantitative 
data with the theory. 

If principle (5) is to hold quite generally, then it must also de- 
scribe the process of the beginning of life from the non-living. Since 
in this process we must obviously have a borderline case, therefore 
we shall have to consider the general expression (5) in which 7 F and 
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H are comparable. In the latter stages of development of the organic 
world the H term will become negligible, and we shall have the re- 
stricted form given by equation (1). This restricted form should not, 
however, be able to account for the origin of life. This throws some 
light upon the difficulty which we encountered in this respect in loc. 
cil. (p. 54). 

In actual applications to borderline cases, we find that friction 
forces play a prominent part. In this case, even in purely physical 
systems, the Hamiltonian principle cannot be applied in its form (2). 
Instead, denoting by f;(dxz;/dt) the frictional forces which are pro- 
portional to the values of dx;/dt , we have (Whittaker, 1927, p. 249): 


x 


dx; ae 
a 6 “;)dt=0. (6) 


J on+sh 


If, as is practically always the case in systems that interest us 
here, the inertial terms are negligible, H depends only on the x;’s, and 
the principle (6) leads to equations of the form: 

dx; oH 
4 + — 0 ; (7) 
dt 0 Xi 


which express the equality of the friction forces and those derived 
from potential H. An estimation of the terms of F,, that depend on 
dx;,/dt indicate, as will be shown elsewhere, that they are in actual 
cases vanishingly small as compared with terms depending on «;’s. 
Hence the generalized Hamiltonian principle which covers both physi- 
cal and biological phenomena will be of the form 


3 ax; 
6 Gé6F —6H — > f,—62,)dt—0. (8) 
x i dt 
That principle leads to the equations 
ax; oH oF 
—— 0 


T (9) 
dt BI05 OD: 


4 


Equation (9) shows that in the general case, 7 0 F/dx; plays for- 
mally the role of a force derived from a potential. If F is small, we 
have the equations (7) of physics. If F is large, we have, for purely 
biological systems, the equations \ 


ax; oF 


=0. (10) 
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The intermediate cases are described by equations (9). 

The sharp distinction between living and nonliving thus vanishes. 
At the same time, however, we neither claim to reduce biology to 
physics, nor the reverse. 
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The expression for the form of organisms in terms of their biologi- 
cal functions, outlined in a previous paper, is carried out somewhat more 
in detail. Making definite assumptions about functions left arbitrary in 
the previous paper, expressions for different shapes of plants are de- 
rived. In continuation of the previous paper, the relation between the 
ee of nervous tissue and the form of an animal is discussed in more 

etal. 

Plants. We have shown previously (Rashevsky, 1943; referred 
to as loc. cit.), how in principle the gross shape of a plant may be de- 
termined by the values of its total mass M , and average rate of metab- 
olism g per unit mass. As a further illustration of the method, but 
without any claim as to immediate applicability of the results, we 
shall, using the same notations as in loc ctt., derive explicit expres- 
sions by making specific assumptions about the functions f(r, 6), 
fo(%,6,M) and f,(qM , 6), which were left undetermined. 

A very rough idea about f may be obtained in the following man- 
ner. A horizontal rigid rod of length 1, fixed at one end and loaded 
at the other end with a load P, will sag at the free end by the amount 
(Geckeler, 1928) 


Lo = 2 B, (1) 
3EJ 
where # is the modulus of elasticity and J is defined by 
J =f ade. (2) 


Here x denotes the coordinate in the direction of the force P, and 
do is the element of surface of the cross section. Hence the dimen- 
sionality of J is [em*]. From equation (1) we find the bending per 
unit length to be 


te 
A LoS Ry (3) 
L 38EJ 


and from the above dimensional considerations, this quantity will be, 
in general, of the form 
69 
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PP 
A % « ——. (4) 
3hr* 
On the other hand, considering that the rod bends under the in- 
fluence of its own weight, we have 


IP ee they? (5) 
Hence 
136 


At <= 
3H r? 


E is likely to increase with 6, and therefore, as a first approximation, 
we may put, assuming H « 6: 


IB 
A Xo « a (6) 


When A 2 exceeds a certain constant determined by the strength of 
the wood, the branch breaks. That constant, in general, increases 
with 6. Let us put it proportional to 6”. Then we have 


l=a 6? 77/8, (7) 


a being a coefficient. This determines f . 

As regards f, , we may consider that the trunk is compressed by a 
force M/x 7,2 per cm? of cross section, and if that quantity exceeds a 
limit, which we again set approximately proportional to 6”, the wood 
breaks under compression. Hence, b being another coefficient, 


M=bD 0? 7,2. (8) 


The function f, may be determined roughly by considering that the 
total metabolic flow in the trunk, which is proportional to gM, is pro- © 
portional to the cross section, and approximately inversely propor- 
tional to 6. Hence, with ¢ as a coefficient, 


Ch ne) COLE Cr 


M = ; =—, is) 
: 6 nN 6 o 
Together with 
M = 6(hr.?2 + nlr?) (10) 
and [loc. cit., equation (27) ] 
qM = kulr, (11) 


equations (7), (8), and (9) determine 5,l,,7,1,7, and n in terms 
of gandM. 
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Solving those equations, we find 


G 1/(p+1) C (p-1)/(p+1) ark? C p/(p+1) 
6= ) ; L=b ae 
bq bq b?q? \ bg ; 


M: q?/2 (PH) azke / @ \P/ (p+) 
as D1/2(p+1) @p/2 (P+) ; =a ; (12) 
q 
see ak? b (6P+5)/ (+1) Mq‘*?+9)/ (p+) 
ia bq?’ hi GS Je eP/ (p+t) 


If » > > 1, then J, , as a function of q, is negative below a cer- 
tain value q, of q , then becomes positive, reaches a maximum, and, re- 
maining positive, decreases to zero, for g = o. A negative I, , physi- 
cally, would probably mean J, = 0, that is, a form of a bush. If the 
actual range of variation of q is in the neighborhood of qi, then J, 
increases with q rather rapidly. For p << 1, l, is first negative, then 
increases to + oo. Here again J, increases very rapidly with q in the 
neighborhood of g,, though for large q, approximately |, « q. As- 
suming that the actual range of variation of q is in the neighborhood 
of g,, we find that the ratio J,/7> increases very rapidly with q. If 
u(q) denotes a rapidly increasing function of q, then, 


lo 

= ee wae (13) 
To VM 

Seen (14) 
Lo qu(q) 

n« Mg, (6<p<7). (15) 
ae aes (16) 


T, @gtl2/M 


For a constant M , small q gives us a small 1/7”, large I/l, , large 7/70 , 
and small ” , hence a shape like in Figure 1a. An increasing q increases 
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1,/%>, and n, but decreases 1/1, and 7r/7.. We have the shape of Fig- 
ure 1b. An increase of M increases n , and to a lesser extent decreases 
l,/7> and r/r,. Hence a.large M, and sufficiently large q give us the 
shape of Figure lc. The density 6 decreases with increasing q , but if 
» >> 1, the variation is small. 

All the above serves, of course, only as an illustration. It indi- 
cates clearly, however, that a theory may be developed on that basis, 
by considering more exact expressions for the different functions in- 
volved. It is important, however, that the expressions (12)-(16) are 
in principle verifiable experimentally and suggest definite new experi- 
mental problems. 

Animals. In loc. cit. we have discussed some possible relations be- 
tween the shape of an animal and the size of its central nervous sys- 
tem. 

We shall conclude this paper by a few remarks in regard to the 
quantitative relations between the values 7,, 72, n,, and ,, intro- 
duced in loc. cit. 

The quantity n is measured approximately by the amount of ner- 
vous tissue necessary to produce and control movements of a given 
complexity. We may approximately put n proportional to the number 
of mechanical degrees of freedom of the lever-propelled metabolizing 
system. Denoting that number by m, we shall have, with A as a co- 
efficient 


n=Am. (17) 


But A itself will depend on at least two factors: one—the muscular 
tonus f to be maintained in each muscle operating on one degree of 
freedom, and two—the amount of special coordination S needed to 
maintain postures which are mechanically unstable. Thus, with B as 
another coefficient 


n=B(t,+S)m. (18) 


For simplicity, let us put t = const. For a simple lever system 
represented in Figure 1a of loc. cit., m =m, is very small and S=0. 
Hence a small n = n,. For a quadruped m = m, is larger, but S = S, 
is still rather small. Hence n,= B(t,. + S,) m, is larger than n,. For 
a biped, m, lies between m, and m,, but S = S, is very large, and 
therefore n, = B(t, + S2)m. > 1,4. 

The quantity m can be calculated in principle relatively simply. 
The calculation of S presents greater difficulties. We may have to de- 
fine properly the degree of mechanical instability of a system, and 
then set S approximately proportional to that quantity. 
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We may now state some of the previous conditions more precisely. 
For instance: 

If v < q/c,andn s Bt, m,, then we have an organism without 
extremities, propelling itself by crawling. 

The case v < g/¢c,n < B ty m, should not occur in nature. 

Ifv > q/c,n = B(t& + Si)ms, we have a quadruped, the length 
of whose legs are determined by equations given in loc. cit. 

Ifv>q/c,n> B(t + S2)m. we have a biped. 

All the above statements are quantitative and can be checked ex- 
perimentally. 

The author is indebted to Dr. H. D. Landahl for reading the paper 
and for valuable criticism. 
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_. The case of a cell in which metabolites are produced and consumed 
within the cell, without leaving the cell, is considered. Conditions are 
given under which there may be an excess production at the center over 
that at the periphery, and this condition is stable with respect to dis- 
placement of the center. The resulting diffusion forces are considered, 
and an expression is developed for the rate of elongation with time. The 
expression is very similar to that obtained by considering the flow of 
metabolites into and out of the cell. 


In previous studies, a theory of cell division has been developed 
on the assumption that there is a net outward dividing force as a re- 
sult of the interchange of metabolites between the cell and its sur- 
roundings (Rashevsky, 1940a; Landahl, 1942). In some cells there 
is very little interchange, while in others division takes place even 
though the various metabolic processes are markedly slowed down 
(Brachet,; 1934). Thus it seems desirable to see if one can modify the 
theoretical mechanism so as to make it essentially independent of the 
interchange of materials with the surrounding medium. 

We shall keep the mechanism as simple as possible and consider 
the effect of a single substance being produced at a rate gq gms cm? 
sec. If we wish to consider quasi-stationary states, it is necessary 
that the substance be consumed in the cell, as we shall neglect the flow 
out of the cell. Further, in order to have any concentration gradients, 
the production or consumption rate must vary with position. Suppose 
that q is constant throughout the cell. This may be the breaking down 
of large molecules into many smaller ones. Then let the rate at which 
the substance resynthesizes or disappears be proportional to the num- 
ber density of a particular catalytic particle, as in N. Rashevsky’s 
polar cell (Rashevsky, 1940b). Let the particles be free to move so 
that their distribution will be given by a Boltzmann distribution. We 
shall neglect the effect of gradients of concentration of the large mole- 
cules on the distribution of the catalyst as well in connection with the 
dividing force. Actually, one might expect several substances to be 
produced and consumed within the cell that also show concentration 
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gradients of the typé discussed here. Not all would have the same 
effect and the resultant would be a weighted average. If the concen- 
tration were higher in the center than at the periphery, the catalytic 
particles would be denser near the periphery, thus tending to main- 
tain the difference. We shall consider this in greater detail. 
Consider a cell which has an axis of symmetry, the z-axis. Let 
2r, be the length along the z-axis, and 27, the width along the z-axis 
perpendicular to the z-axis. We shall treat the problem by using the 
approximation method (Rashevsky, 1940). Divide the cell as in Fig- 
ure 1, so that V,; = V; where V; is the volume of section 2, and V is 


—— 


FIGURE 1 


the volume of the cell. If c; is the average concentration and 7, is the 
average number density of the catalyst in section 7, we have 


N, Vy + Ne Ve + Ng Ve= NV, (1) 


and a similar expression for the c;. In equation (1), ” is defined by 
the equation as an average value of the number density of the catalyst 
for the cell as a whole. 
If akT/N is the constant of proportionality between the force on 
a catalytic particle and the negative of the concentration gradient, we 
have 
Ti pet ed, (2) 
to determine the distribution of the catalytic particles. R is the gas 


constant per mole, T is the absolute temperature, and N is Avogadros’ 
number. 


We next determine the conditions for the average concentration 
to be higher in section two than that in one and three. From the ma- 


terial balance conditions, we have, denoting by d; the distance from 
the center of the cell to the center of section i, 
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dc, AD 

Fy eran ean Se 
dc. D AD 
ee ies V, (ce =-"e7) Ge (Coe —3Cz) & (4) 


and an expression for section three equivalent to (3) with subscript 
1 replaced by subscript 8. A is the area of the plane separating sec- 
tion two from three or one. The constant a determines the rate of re- 
moval of the substance by the catalyst. First we shall let d, = d,, 
that is, assume that the concentrations are symmetric with respect to 
the z-axis, thus c, = c; , and similarly n, = n;. Introducing equations 
(1) and (2) in (4), we have 


dc. 2AD anV 
(c2 — Gi) 


= — (5) 
dt d, V. PANE: e2(C2-€1) AP V. 


If de./dt = 0, and (c — ¢c,) = 0, equation (5) is satisfied, since 
q = an from the definition of n in equation (1). We wish to know 
under what conditions dc./dt > 0, for ce — ¢ > 0, cc. — ¢, being 
small. This condition is readily seen to be 


aV,V.dq>AVD, (6) 


that is, the force on the particle per unit gradient of concentration, 
or the production rate q, or the size of the cell must be large enough, 
or the internal diffusion coefficient D must be small enough, in order 
that a small increase in c, will cause a further increase. Similarly, 
we find that a small decrease in c, will cause a further decrease if in- 
equality (6) holds. Thus (6) is the condition for the equilibrium 
C; = C, = ¢; to be unstable. 

If a(c. — ¢,) is sufficiently large, so that we may neglect the term 
containing e*) in the denominator, we may solve for (cz, — ¢,) in 
equation (5) for the equilibrium value by setting dc./dt = 0, obtain- 
mg, 10r ¢e.— > 0, 

Ved. q 
2AD- 


If V, © 4V, ~ 1 J, this is a(¢, — ¢,) = 2aV,V2d,q/A DV > 2 from 
(6). Then also, (7) may be made exact if the right side is multiplied 
by tanh 4 a(c. — ¢:) > .76. If the left member of (6) exceeds the 
right by a factor of four, the correction factor differs from unity by 


but a few per cent. 


(7) 


Cola Ci 
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Suppose next that the point of maximum in section 2 of Figure 1 
is displaced slightly to the right, so that we may treat this effect ap- 
proximately by introducing d, > ds. We then wish to know whether 
the average concentration in section 3 will tend to rise or fall. If it 
rises, the maximum would tend to move to the right, and if this con- 
tinued, we would have a polar cell (Rashevsky, 1940b). On the other 
hand, if ¢; falls, the maximum would be shifted back and the sym- 
metric configuration can be considered stable. 

From (3) and (2), we have for the stationary state 


AD No 
log —=0, (8) 
1 Vio NM 


q— a, — 


and a similar expression for section 3. Subtracting (8) from the cor- 
responding equation with subscripts 1 replaced by 3, introducing 
6= (d, — ds) /d, < <1 , and'expanding loge 2,/", for 1, — 2%: —— te, 
we have on solving for n; — n,, 


AD 6n; log Ns3/Nez 


— t= : 
+ : an;oad,Vi,—- AD 2 
Since we are taking 6 > 0, c > c¢;, thus n; > m, and a ViVed.q > 
ADV from (6), then, > .,as am; > @, V > Vez. But % > ine 
plies that cs; < ¢, or c; falls below c, , thus below its original value for 
=0. Thus the central maximum is a stable equilibrium if (6) holds. 
Heretofore we have considered the stability with respect to the 
z-axis. In this case we found (¢ — c,), to be given by, equation (7). 
In a similar manner we could treat the stability with respect to x and 
obtain (c; — ¢,)¢. However, in this case ¢,, would not necessarily be 
the same as ¢., , whereas physically this must be so. 
If we divide the cell, as in Figure 2, into a cylindrical central 
region 2, of volume V,, a region 1 of volume V, surrounding region 


FIGURE 2 
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2 coaxially, and two regions comprising the “ends” of the cell and of 
volume V,; together, we may write an equation formally the same as 
equation (1), and two independent equations of the form (2). We 
assume a center of symmetry. From these we may solve for »; in 
terms of € and ¢ where 


E=al(C.—c,), C=a(c. — ¢s). (10) 


The ¢;’s are concentrations corresponding to the new divisions. For 
simplicity, set V; = V, = V; = V/3. We may also write three equa- 
tions similar to (8) and (4). Note that A, # A; and q= an. Sub- 
tracting the first of these equations from the second, then the third 
from the second we have 


dé 
— = 8ag(e—1)/(1 + ef + et) — 2K 8 — Kit, (11) 
0 

dc 

Gp eae — V/A + + et) Shee Kec; (12) 


where 
K,=3A;D/d,V. (13) 


As K,~© K,, and if agdV > 3KdV =9 A D, which is essentially 
the same as (6), then dé/dt > 0 and df/dt > 0 if € ~ ¢ are small posi- 
tive values. Thus & and ¢ continue to increase. If K, is nearly equal 
to K,, and K is a value between K, and K,, then for e& > > 1, 
&=C=aq/2K,isa solution for the stationary state. Now let 0 —¢ 
— &€ << &. We shall obtain a first approximation to © , expanding 
e® and retaining the first degree terms. Letting d&/dt = df/dt = 0, 
and eliminating ¢ by ¢= & + ©, we have approximately 


ee (K, — Ks)é 
Bk ke dg Ke 


Introducing & = aqg/2K , ag > > 2K/3, and 6K > > K, — K;, equa- 
tion (14) becomes 


(14) 


hon 
6@= Ss (15) 
6K 


We consider next the influence of the concentration gradients on 
the shape of the cell and follow the approximation treatment devel- 
oped by N. Rashevsky (1940a, Ch. iii). We assume that there are 
structures within the cell upon which the diffusion forces may act 
whose relative volume is », and that the force per unit volume is given 
by ; 
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F=-- grad c. (16) 


The relative elongation rate (dr,/dt)/1r1, 2r, being the “length”, is 
given by integrating over the surface and volume forces. The first 
is represented by the surface tension y and contributes an amount 
(Rashevsky, 1940a, Ch. iii, eq. 24) 


“a y(r — 2) /2%1 To (17) 


where 27, is the width of the cell and 7 is its viscosity. The contribu- 
tion from the volume forces is (cf. Rashevsky, 1940a, Ch. iti, eq. 9) 


—~3RT u(cs— ¢:)/2Mn, (18) 


where a factor of two is introduced here, because c; and c, are taken 
half way out from the center of the cell instead of at the “ends”. 
Adding the expressions (17) and (18), introducing (15), (18) 


and: (65: -¢;)) = —1O/a.; with Aso 17, 74/2 Ag tie 
A, —ar? ,d,—72/2, ds=7,/2, and setting 
e= (1, — 7%) /%; (19) 
we have 
de 
ae ee (20) 
with 
SES 
re fe, , (21) 
1 3y 
B=—-—, 
47, 4n -) 
1 ay 
C= —8K'’ +—), 
ae e (23) 
,__ 8RT pro 
4M no : (24) 


3 
If we set a = 5 Rk T V,/M, so that the force on each catalytic particle 


of volume V, is given by equation (1) of N. Rashevsky (1940a, Ch. 
iii), we have 


ast. RT ur 
27V,N- 


v 


(25) 
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Note that K’ is independent of q. If one assumes that e® is large com- 
pared to one, K’ becomes proportional to q, but the above may not 
hold. 

Equations (20) to (23) are very similar to the equations result- 
ing from the case of an outflowing metabolite (Landahl, 1942, eq. 6— 
9). The derived parameter K’ and y/y are of the same order of mag- 
nitude as in the previous case. Note that C < 0, so that the slope 
curve (Landahl, 1942, Fig. 1) would have its maximum displaced 
slightly to the left of the middle, whereas the measured values have a 
maximum just to the right. 

If we take K’= 7 X 10 (Landahl, 1942), w= 10", andy» =}, 
we obtain for the radius of the catalytic particles, 7, = 1.6 X 10-> em. 
Then, if M = 10? ,a=2 X 10’, so that the inequality (6) becomes 


D <20q. (26) 


The rate of accumulation of lactic acid in developing Arbacia eggs is 
of the order of 3 X 10°° gm cm sec", (Hutchens et al, 1942) and the 
corresponding q must be at least as large. Thus if we consider lactic 
acid as being formed from large glycogen molecules, but also being 
resynthesized by the action of a catalyst, we have a special case of the 
mechanism assumed initially. Taking q to be somewhat greater than 
the rate of accumulation, g ~ 10, then the diffusion coefficient for 
lactic acid within the cell should be of the order of 2 X 10-*. But the 
value of the diffusion coefficient in water is 1.4 X 10-° cm? sec? . Thus 
the ratio of the diffusion coefficient inside the cell to that outside must 
be less than one-seventh, or a correspondingly larger fraction if q is 
larger. But since these molecules are polar, the diffusion coefficient 
within the cell would be expected to be very much smaller than that 
in water, so that the ratio would be much less than one-seventh. Thus 
inequality (6) is satisfied for this special case. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of The University of Chicago. 
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